Fluctuation induced first-order phase transitions in a dipolar Ising ferromagnet slab 
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We investigate the competition between the dipolar and the exchange interaction in a ferromag- 
netic slab with finite thickness and finite width. From an analytical approximate expression for the 
Ginzburg-Landau effective Hamiltonian, it is shown that, within a self-consistent Hartree approach, 
a stable modulated configuration arises. We study the transition between the disordered phase and 
two kinds of modulated configurations, namely, striped and bubble phases. Such transitions are 
of the first-order kind and the striped phase is shown to have lower energy and a higher spinodal 
limit than the bubble one. It is also observed that striped configurations corresponding to different 
modulation directions have different energies. The most stable are the ones in which the modula- 
tion vanishes along the unlimited direction, which is a prime effect of the slab's geometry together 
with the competition between the two distinct types of interaction. An application of this model 
to the domain structure of MnAs thin films grown over GaAs substrates is discussed and general 
qualitative properties are outlined and predicted, like the number of domains and the mean value 
of the modulation as functions of temperature. 



I. INTRODUCTION 

Magnetic phase transitions in materials with finite spa- 
cial dimensions is still a subject with many aspects to be 
understood and to be investigated theoretically. In these 
systems, one finds the competition between a strong, 
short-range interaction (exchange) and another weak, 
long-range one (dipolar), from which a modulated sta- 
ble configuration is expected to outcome Q. Due to the 
variety of possible modulated patterns, like striped, bub- 
ble and intermediate shapes, many different complex do- 
main structures are Hkely to be seen. This phenomenon 
is observed not only in magnetic materials 0, 1^, but 
also in other systems characterized by the same kind 
of competition between an organizing local interaction 
and a frustrating long-range interaction for example, 
spontaneous modulation of mesoscopic phases is found 
in biological systems, amphifilic solutions |4| , Langmuir 
monolayers [3 and block copolymers ^QJ. 

There are, indeed, several analytical and numerical 
studies in the literature investigating size effects on the 
critical behaviour of magnetic systems 0,&0|- However, 
the majority of them deals with periodic boundary con- 
ditions, and not free boundaries, which comes to be the 
case for many materials. Even the works related to this 
latter kind of systems ^3, ^| do not take into ac- 
count the weak, long-range (dipolar) interaction, which 
changes drastically the underlying physics. As shown by 
Garel and Doniach the inclusion of the dipolar in- 
teraction in the Ginzburg-Landau effective Hamiltonian 
of a magnetic slab with infinite width leads to a min- 
imum in the Fourier space characterized by a non-zero 
wave vector. This is responsible not only for instabil- 
ity towards the spacial homogeneous phase but also for 
the existence of a large volume, in the Fourier space, for 
fiuctuations of the order parameter to take place and in- 
duce a first-order transition (Brazovskii transition [T^l. 
Therefore, to achieve a more complete understanding of 



finite magnetic systems, it is necessary not only to con- 
sider the finiteness of the them but also the fiuctuations 
of the order parameter. 

To explain some properties of many real materials, size 
effects are in fact necessary. For example, concerning 
non-magnetic systems, Huhn and Dohm have proposed 
that size effects are responsible for the temperature shift 
of the specific heat maximum in confined He'* In 
what concerns magnetic systems, a material that has 
been deeply experimentally investigated in recent years 
and in which size effects may play a significant role is 
MnAs thin films grown on GaAs [lj|. The reason why 
this heterostructure is caUing so much attention is due 
not only to its academical appeal but also to its possi- 
ble application as a spintronic device In contrast 
with bulk MnAs, which presents an abrupt transition 
from the low temperature hexagonal (ferromagnetic) a 
phase to the hi gh t emperature orthorhombic (paramag- 
netic) /? phase [13, the MnAs: GaAs films show a wide 
region of coexistence between a and (3 phases from ap- 
proximately 0°C to 50°C [Hill Ea El. In this 
region, periodic stripes subdivided in ferromagnetic and 
paramagnetic terraces arise. X-ray diffraction experi- 
ments i2(3] and microscopy measurements have shown 
that, while the temperature varies, the width of the ferro 
and paramagnetic terraces change, but the stripes re- 
main with the same periodic width. These experiments 
have also brought out the terraces morphology, showing 
the complex phases inside the ferromagnetic terraces. As 
their width is of the same order of magnitude as their 
thickness, both spacial limitations are important to un- 
derstand their internal domain structure. 

Here is an outline of the article: in Section 2, we con- 
struct an expression for the Ginzburg-Landau effective 
Hamiltonian of a dipolar Ising ferromagnetic slab with 
finite width and finite thickness, considering Dirichlet 
boundary conditions (i.e., vanishing of the order param- 
eter at the walls). We show that a modulated phase 
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arises as the ordered one, and is represented by a dotted 
semi-elUpsis in the Fourier space as a result of frustra- 
tion and Dirichlet boundary conditions. In Section 3, we 
apply a self-consistent Hartree calculation to take into 
account the fluctuations of the order parameter around 
the region of minimum energy. Generalizing the original 
method developed by Brazovskii to the case of this 
flnite system, we calculate the free energy proflles for 
two different types of modulation: striped and bubble 
phases. We show that the striped phases are more stable 
than the bubble ones and also that the energy degeneracy 
concerning the region of minimum energy is broken along 
the flnite direction. In Section 4, we discuss the applica- 
tion of the model to the real case of MnAs:GaAs films. 
Although in such systems the magnetization is rather 
vectorial than Ising-type, general qualitative properties 
due to the nature of the interactions and to the slab's 
geometry can be obtained. Section 5 is devoted to the 
final remarks and followed by an appendix where details 
of some calculations are explicitly derived. 



Tc is the Ising ferromagnetic critical temperature, a is 
the lattice parameter and m{r), the scalar order param- 
eter, is the coarse-grained spin in the position f— {x,y). 
We are denoting, through all this article, the integrals 
over the region limited by the plane of the slab as: 

f 2 f°° fd 
d r ^ / dxdy . 

J J -oo Jo 

For the sake of simplicity, we consider a cubic lattice. 
In the long wavelength limit, the actual crystallographic 
structure will not change significantly the basic physical 
properties of the system. In this limit, we can calculate 
the dipolar contribution to the total energy fdip using 
just Maxwell equations. 

To obtain fdip, we express the magnetization M{r) at 
the position r in terms of its Fourier components as 



II. GINZBURG-LANDAU FOR THE FINITE 
SLAB 



M(f) 



E 

n>0,qy 



sin (—J e-'^yy , (3) 



We consider a slab with thickness D {z axis), width 
d (x axis) and no limitation along the y axis; the mag- 
netization is supposed to point only to the z direction 
(Ising model) and to depend upon x and y only (uniform 
along z). This is, indeed, a very simplified model of a fer- 
romagnetic stripe on the MnAs:GaAs coexistence region, 
but we will postpone this discussion until the last section. 
Our purpose in this section is to obtain a two-dimensional 
Ginzburg-Landau to describe the system; first, consider 
the well-known mean field expansion of the free energy 
due to the exchange interaction ^22] : 



(1) 



where 
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where g is the gyromagnetic factor, /zb is the Bohr mag- 
neton and: 
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i{r) sin ^- 



d I 



e'lyVdxdy . (4) 



In expression (EJ, the sine term appears as a conse- 
quence of the boundary condition that the magnetization 
vanishes at the edges of the slab (Dirichlet boundary con- 
ditions). In what concerns MnAs:GaAs thin films, this 
condition approximates the fact that, in the coexistence 
region, the ferromagnetic terraces are succeeded by para- 
magnetic ones. 

As we show in Appendix A, the magnetostatic energy 
of the arbitrary configuration IpJ can be straightforward 
calculated from Maxwell equations, yielding: 



du 



(l _ e-pV«'+(««'^)') [l + (-1)™'+! cos(u) 



+ (qyd)^ (u^ — n^TT^) (u^ — n'^Tr^) 



(5) 



where the sums are over n and n' with same parity and properties of the slab, we look for an analytical approx- 
where we have denoted the slab's aspect ratio asp = d/D. imation for Disregarding the cross terms, which are 
As we wish a simple model to describe the main physical usually negligible, we have that the sum of the direct 
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terms is: 

fMp = [— ) 2^ TT — m„,,^m„,_,„ , (6) 

n,qy 

where we denoted the wave vector modulus by: 

1=\l^+1y- (7) 

The accuracy of the approximation ij^J depends on the 
values of the parameters involved and on the pair {n, qy) 
considered. In the experimental case of interest, namely, 



the MnAs:GaAs films in the neighbourhood of the phase 
transition between the ordered and disordered phases, 
this approximation implies in errors less than 20% as long 
as p > 0.5. 

The approximate analytical expression obtained is very 
similar to the expression deduced by Garel and Doniach 
for a slab with infinite width The only difference 
is that, in the present case, the wave vector component 
along the x direction is discrete due to the slab's finite 
width. It is clear that in the limit d ^ oo we recover the 
same expression. 

Using equation © , we obtain the following expression 
for the total free energy density 
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The quartic term is the same as those obtained in pre- 
vious works about finite systems with Dirichlet boundary 
conditions (see, for instance, [1(1 H^l. The effects of the 
two interactions mentioned before are evident from {HJ: 
while the q^ term, generated by the exchange energy, 
favours q — (non-modulated) configurations, the last 
term, generated by the dipolar energy, favours q ^ oo 
configurations. The total energy reaches its minimum 
value when the wave vector modulus is given by (as long 
as qoD » 1): 



1 ^ 16^.9^1 



1/3 



(11) 



which means that the most relevant thermodynamic 
states are characterized by a non-zero modulation {qo ^ 
0). Expression lfTT|l is the same as the one obtained by 
Garel and Doniach in for the minimum energy of the 
infinite slab; however, in their case, the phase space of 
low energy excitations is described a circle in the Fourier 
space, whereas in the present situation it is described by 
a dotted semi-ellipsis in the {n,qy) space, as it is shown 
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Figure 1: Region on the momentum (Fourier) space corre- 
sponding to the minimum of the Ginzburg-Landau for the 
case of a slab with (a) infinite width and (b) finite width. 



in figure ^ 

Expanding JSJ around its minimum, it is straightfor- 
ward to obtain the respective partition function: 
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Z = J Vmexp{-n[m]) , (12) 

where the Ginzburg-Landau effective Hamiltonian Ti. is 
given by: 



1 



+ - / d^rm'^ir) , (13) 



n[m] ^ J d^rcf r'm{r)GQ ^(r, r)m{'r ) Fourier series as: 
I 



and the correlation function Go is written in terms of its 




The parameters u, tq and c appearing in the above 
expressions can be written in terms of the microscopic 
parameters of the systems as: 



[T -T*)D 



8aT 



(15) 



where we defined the shifted critical temperature: 



T* = Tr 
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1/3- 



(16) 



Therefore, considering the expression for the correla- 
tion function, Ijl4|l . it is expected some similarity between 
this system and the Brazovskii's model The question 
is if the reduction in the momentum (Fourier) space, due 
to the discreteness of the x components {qx = n-K /d), is 
able to substantially change the picture, since Garel and 
Doniach have shown that, for the case of an infinite slab, 
in which qx is continuous, a fluctuation induced flrst- 
order phase transition does occur between the ordered 
(modulated) and the disordered phase 0- 



III. HARTREE CALCULATION AND 
BRAZOVSKII'S PROCEDURE 

As we detailed in the previous section, there is a de- 
generate region at the Fourier space, expressed by the 
non-zero wave vector modulus go, corresponding to the 



minimum of the Ginzburg-Landau effective Hamiltonian. 
Hence, there is a large space for fluctuations of the or- 
der parameter to take place, and a mean-field approach 
to calculate the partition function lfT2ll - and its corre- 
sponding thermodynamical properties - is not satisfac- 
tory. To deal with that, we generalize the procedure 
adopted by Brazovskii to our finite system. Such 
procedure is based on the Hartree self-consistent method, 
that consists on replacing the quartic term by an effective 
quadratic one |2^ : 



J Jd\m\r)-^Y J d\ {m\r)) m\r) . (17) 
Then, using the identity 

(m^{r)) ~ G{r^r) + {m(r))'^ , 

and substituting in II13|I . it is clear that a self-consistent 
equation is obtained for the correlation function: 

G-^ry) = Go\r,r') + 3uG{r,f)S{f-r') + 

+3u{m{r)f Sif-f*) , (18) 

where the Dirac delta function is to be understood as 
belonging to the intervals [0, d] {x axis) and [— oo, co] {y 
axis), and not to [—00,00] and [—00,00], as it is usually 
assumed. 

Following Brazovskii, the dominant contributions for 
the correlation function G(r, r') come from the di- 
agonal Fourier components G^n.qy),{n.-qy) such that 

^Jr}^TT^J(P~i~q^ — qo. Therefore, it is useful to write 

the self-consistent equation on the Fourier space 
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where 



and their mean values can be evaluated to yield 



(20) 



The summation in Ijl9ll is over the region compre- 
hended by the dotted semi-ellipsis shell whose thickness A 
is such that A ^ go- In this region, the diagonal Fourier 
components Gn,g„ can be expanded as: 
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Here A^(x) denotes the integer closest to x and smaller 
(21) than x. Therefore, equation l(T9|l can be written as: 
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where: 
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It is clear that the self-consistent equation for the 
renormalized parameter r , equation Ij23|l , depends on the 
system phase through {m{r)). However, for any configu- 
ration {m{r)), analogously to fl^, the equation does not 
allow the solution r = 0, what implies that the transition 
from the disordered to this ordered phase (characterized 



by the mean value {mif))) is not second-order. There- 
fore, a first-order transition is expected, as well as the 
raise of metastable states (and the respective spinodal 
stability limits). 

To further investigate how the system achieves all the 
possible distinct modulated states, it is useful to calcu- 
late the free energy difference between these configura- 
tions and the disordered one. Hence, we use the same 
procedure as Brazovskii and consider that the conjugate 
field h grows from zero in the disordered phase to a max- 
imum value and then goes again to zero in the ordered 
stable phase, characterized by an amplitude A ^ '■ 
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(25) 



where 
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SF 



is the conjugate field in the Fourier space. A lengthy but straightforward calculation gives: 



m.m ,p,p 



m^n,p 



3u 3w 



(26) 



We follow Garel and Doniach j3| and study two differ- 
ent modulated configurations: the striped and the bubble 
phases. 

Striped phases 

The name stripes may be somehow misleading, as gen- 
uine stripes cannot be formed due to the boundary con- 
ditions, unless they He only along the y direction. In fact, 
these phases refer to the simplest modulated configura- 
tion that can arise in the finite system: 
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A5n,no \^q,qay + ^q,~qoy 



(m (r)) = 2ylsm \—^ j cos (go^y) 



(27) 
(28) 




(a) 



(b) 



Figure 2: Contour plot of the order parameter along the slab's 
plane concerning the simplest modulated phases that arise in 
the case of (a) an infinite slab and (b) a finite slab. 



where 



qo ■ 



Therefore, the use of the name stripes is to make con- 
tact with the phases of the infinite slab rather than 
to describe exactly the geometrical pattern. Figure [2 
compares the general picture of the simplest modulated 
phases of the infinite and of the finite slab. 

Substituting equation Ij28|l in expressions II23|I and 
II26II . we obtain the following expressions for the self- 
consistent and the state equations: 



r = ro + 



-uA' 



(30) 



This equation is the same as the one obtained for the 
case of an infinite slab, and was considered by Brazovskii 
in his original work. It implies that these striped phases 
can arise as metastable states below the spinodal stability 
limit: 

^spinodal ~ 1.89 (wF^q ) 

To obtain the free energy difference between these 
modulations and the disordered phase, we take (m (r)) — 
in we get: 



h = rA uA-^ 
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(29) 



Imposing that the conjugate field vanishes in the or- 
dered phase (denoted by r^), we get: 



r = ro + 



(31) 



Using II25|I ■ the free energy difference is 



Figure 3: General behaviour of the free energy difference be- 
tween the striped phase (no = 1) and the disordered phase 
as a function of the ratio qod/n. The temperature in (a) is 
greater than the temperature in (b). The shaded regions in- 
dicate that the ordered state is not metastable (temperature 
above the spinodal) and therefore it does not make sense to 
define a free energy difference. 



Figure 4: General behaviour of the free energy difference be- 
tween the ordered and the disordered states concerning the 
striped phase with no = 1 (full line) and the striped phase 
with no = 2 (dashed line). The shaded regions indicate that 
the ordered state is not metastable (temperature above the 
spinodal) and therefore it does not make sense to define a free 
energy difference. 



2{uTt 



? 9 

Pa P I — 



(32) 



where we defined the auxiliary variables 



Pi 



Ti 



("r„o) 
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The new feature that appears as a consequence of the 
finiteness of the system is the dependence of the factor 
r„(, with respect to the slab's width d and to the modu- 
lation label no, which indicates what point of the semi- 
ellipsis is taken to modulate the system. It is clear that 
a modulation labeled by no can arise only if: 



> no 



(34) 



otherwise the semi-ellipsis does not comprehend this spe- 
cific point. Such label can be interpreted as the number 
of "spread" domains along the x direction, since there 
are no sharp walls and the magnetization changes sign 
continuously from one domain to another. 

From the definition of r„Q, equation Ij24|l . it is clear 
that when the ratio qod/n is an integer the factor di- 
verges. However, this does not mean that the free energy 
difference diverges, since the spinodal limit also depends 
upon r„Q . A plot of this energy as a function of the ra- 
tio qod/ir for a given temperature is shown in figure El 
Firstly, it is clear that there are barriers in the energy 
profile whenever: 



is the break of the semi-ellipsis degeneracy. In a mean- 
field approach, all the different modulations that can ex- 
ist for a certain slab's width d would have the same en- 
ergy. However, as shown in figure 0) when we take into 
account the fiuctuations, each modulation no assumes 
distinct energy values, implying that the high degener- 
acy of the minimum energy is broken. Moreover, as the 
slab's width increases {d oo), the energies get closer 
again, what agrees with the result for the infinite slab, 
where there is no degeneracy break. 

A deeper analysis reveals that, in fact, only the degen- 
eracy along the qx direction is broken, and not the other 
along the qy direction. This is only a refiection of the 
translational invariance break along the x direction, due 
to the existence of edges. 

It is important to analyze carefully the energy barriers 
that appear when condition l(35ll is met, because in such 
situation (and only in it) genuine stripes, characterized 
by modulation only along the x direction, can appear. 
This configuration is given by 



(36) 
(37) 



and does not obey to the same Hartree or state equations 
of the configurations previously considered. Indeed, a 
direct substitution of l|H7ll in (EHJ and ll^ imply that: 



ro 



(35) 



We also note that, as the temperature decreases, the 
barriers heights become approximately uniform. Another 
important consequence of the r„Q dependence upon no 



h = rA uA^ . 
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Hence, the self-consistent equation is given by: 



rA' 



= ro 



(38) 



(39) 
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Figure 5: General behaviour of the free energy difference be- 
tween the ordered and disordered states concerning the "gen- 
eral" striped phase (no = 1 - full line) and the genuine striped 
phase with no modulation along the y direction (no = 1 - 
dashed line) as a function of the ratio qod/n. The shaded re- 
gions indicate that the ordered state is not metastable (tem- 
perature above the spinodal) and therefore it does not make 
sense to define a free energy difference. 



direction is shown in figure It is clear that the lat- 
ter is always more stable than the former; however, it is 
important to bear in mind that the genuine stripes can 
only appear when condition is fulfilled. The spinodal 

limit for them is given approximately by —1.5 (uTno)'^^^, 
what means that these configurations always appear be- 
fore the "general stripes". Therefore, we can say that 
whenever the geometric condition is met and the 
system can be divided in non-modulated configurations 
along the y direction, it will do. 



Bubble phases 

In an infinite slab, where the wave vector components 
Qx and qy are continuous, an hexagonal bubble phase is 
described by the configuration 



and the free energy difference by: 



AFv = — ^ — — 
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(40) 



A plot comparing the energy profiles of the "general" 
stripes and the stripes with no modulation along the y 
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and 



go . 



(41) 



It is clear that such condition can no longer be satisfied 
by the finite slab, due to the boundary conditions. Than, 
to study other phases than the simplest "striped" ones, 
we consider a phase that resembles some aspects of the 
bubbles in the infinite slab, as shown in figure 



(m(r)) = Asin (^^^^^^^^ cos((7oyj/) + ^sin 



(no + 1) TTX 



cos{qiyy) , 
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where: 



d2 



(no + ly ^TT^ 



qly = qo . (43) 



Once more, we adopt the name bubbles to keep the 
correspondence to the case of the infinite slab, and not 
to describe the actual geometric pattern. We note that 
this configuration can take place as long as: 



qod 



> no + 1 



(44) 



Substituting expression II42|I in Il2;^|l and Ij26|l yields the 
following Hartree and state equations: 



Since the conjugate field vanishes in the ordered bubble 
state, we obtain, for the self-consistent equation: 



ZTA =ro + 



The spinodal is approximately —2.51 {uT„ 



(46) 



x2/3 



what 

means that this configuration appears after the striped 
phases. In a mean-field calculation, they (and any other 
modulation characterized by go) would appear simultane- 
ously for both the infinite and the finite slab. Therefore, 
the fiuctuations break also this degeneracy, but this is 
not an effect due to the finiteness of the system, since it 
occurs also for the infinite slab. 

Using Ij45ll and Ij46|l , the energy difference is calculated 



rA uA^ . 
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(a) (b) 

Figure 6: Contour plot of the order parameter along the slab's 
plane concerning the "bubble" phases that arise in the case of 
(a) an infinite slab and (b) a finite slab. 




Figure 7: General behaviour of the free energy difference 
between the ordered and disordered states concerning the 
striped (full line) and the bubble phase (dashed line), both 
with no = 1, as a function of the ratio qod/ir. The shaded 
regions indicate that the ordered state is not metastable (tem- 
perature above the spinodal) and therefore it does not make 
sense to define a free energy difference. 



As it is the case for the striped phases, energy barriers 
are observed when condition (i;-i5|l is met and, for differ- 
ent values of the label no, different values of energy are 
observed. Figured compares the energy difference of this 
bubble phase with the one referring to the striped phase; 
we note that, in general, the former is greater than the 
latter, what means that the simplest modulation is more 
stable. Moreover, since the bubble's spinodal is lower 
than the stripe's spinodal, we expect that the magnetic 
domain configuration will never be divided into bubbles. 

This picture can change in the presence of an exter- 
nal magnetic field along the z direction. As showed by 
Garel and Donicah |3| in the case of the infinite slab, in 
a mean-field approach, the magnetic field can favour the 
formation of bubbles instead of stripes for certain tem- 
perature ranges. We expect that, in the present case of 
the finite slab in the Hartree self-consistent approach, a 
similar phenomenon can occur. However, since this is not 
the scope of this work, we do not investigate further such 
subject . 

More complex patterns built up from other combina- 
tions of the semi-eUipsis points are also possible; however, 
the calculations involved become more difficult. From the 



previous analysis, we expect that the simplest modula- 
tion will be the most stable one, as it is the case for the 
infinite slab. The main difference is that, in the latter 
case, the spinodal of more complex phases are greater, 
and not lower, than the spinodal of the simplest modu- 
lation. 



IV. APPLICATIONS TO MnAs : GaAs FILMS 

The aspects presented in the last section are particu- 
larly interesting on systems in which the slab's width d 
can be varied. In fact, this is the case for MnAs thin films 
grown over GaAs substrates, where it is observed the for- 
mation of ferromagnetic terraces whose widths depends 
almost linearly upon the temperature pS]: 

d(r) = 600- 12T, (48) 

where d is given in nanometers and T in Celsius degrees. 
This is valid in the region where the ferromagnetic ter- 
races coexist with the paramagnetic stripes, from °C to 
50 °C. In this section, we intend to discuss the domain 
structures inside the ferromagnetic terraces, considering 
them as finite slabs, and compare to experimental results. 

It is important to notice that, in the MnAs:GaAs sys- 
tem, the spins responsible for the magnetism are not 
scalar (Ising-like) , but vector (due to the crystalline field, 
they would be better described by a xz model, and an 
approach following the fines of pj| would be necessary) . 
Besides, the film thickness is larger than 100 nm, what 
means that three dimensional domains could be formed. 
Nonetheless, as we are concerned with the general pic- 
ture of the problem, we believe that this simple model 
proposed can outline some general properties due to the 
nature of the competing interactions (strong short-range 
versus weak long-range) and to the geometry involved 
(Dirichlet boundary conditions in Cartesian coordinates). 
However, the specific features of the domains that would 
be formed can be much more complex, as we showed pre- 
viously for the case of one-dimensional Neel walls . 

First of all, we need to estimate the order of magnitude 
of the parameters. We do not intend to obtain an exact 
quantitative description, but rather some qualitative in- 
sights about the domain structure of each ferromagnetic 
terrace. Therefore, based on the experimental studies re- 
garding MnAs:GaAs thin films [11 El ElES El, we 
take a = 5A, (7 = 3, = 130 nm and Tc = 32 meV. 

Substituting these values in the equations deduced in 
the previous sections, we can study the behaviour of 
striped and bubble phases inside the ferromagnetic ter- 
races. Figure IHI shows the energy difference between the 
striped phases and the paramagnetic (disordered) phase. 
We note that the temperature is always below the spin- 
odal limit and that there are local energy minima refer- 
ring to the configurations in which there is no modulation 
along the y direction. Such configurations comprehend 
structures from 1 (the last energy minimum) to 10 (the 
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Figure 8: Free energy difference between the striped and the 
disordered phases for MnAs:GaAs films as a function of tem- 
perature (in Celsius degrees). The sharp minima refer to the 
configurations for which there is no modulation along the y 
direction and n domains along the x direction, where n goes 
from 10 (first minimum) to 1 (last minimum). Contour plots 
for some of these configurations are presented, as well as a 
sketch of this profile in the case of a real system. 

first energy minimum) "spread" domains lying along the 
X direction. Note that, for a real system, in which there 
are impurities and the various ferromagnetic terraces do 
not have exactly the same width at a given temperature, 
these local minima would not be so sharp and the free 
energy profile would be continuous, as sketched in the 
figure. 

We also see that, when T w 45°C {d k, 550 nm), the 
semi-ellipsis at the Fourier space representing the mini- 
mum energy does not comprehend any positive and in- 
teger n, what means that the modulated state cannot 
arise anymore. Above this temperature, it is likely that 
the interaction between neighbour ferromagnetic terraces 
will play an important role to determine the new stable 
configurations. 

Another aspect that is not represented expHcitly in 
the figure is that the break of degeneracy between states 
corresponding to different number of domains n is very 
weak, since the energy scales involved are experimentally 
unnoticeable. This is due not only to the large value of 
the slab's thickness D but also due to the fact that the 
temperature is far below the spinodal limit. Therefore, in 
the beginning, when the temperature is just above 0°C, 
all configurations with n < 10 and non-zero qy (such that 
the wave vector modulus is go) are degenerate. As the 
temperature increases, the system meets the first local 
free energy minimum, corresponding to n — 10. Hence, 
at this temperature, the terrace will be divided in 10 
domains along the x direction and no modulation along 
the y direction. In the sequence, all configurations with 
n < 9 and qy ^ Q are again degenerate; however, since 
the system was previously in a 10 domain configuration, 
it will be energetically favourable to "destroy" just one 
domain along the x direction before it meets the local 
minimum corresponding to n = 9. Therefore, we expect 
that each local free energy minimum corresponds to a 
change in the number of domains along the x direction. 




10 20 30 40 50 



Figure 9: Mean value of the MnAs:GaAs films modulation 
along the y direction, in units of 10^ A , as a function of the 
temperature in Celisus degrees. A Gaussian distribution is 
considered for the widths of the sample's terraces. 



and that, for these specific temperatures (where the free 
energy has these local "traps") , the configurations will 
be non-modulated along the y direction. 

An experimental measure that can show the occurrence 
of these local free energy minima is the mean value of qy, 
the y-component of the modulation, that can be obtained 
by x-ray scattering. As we have already pointed out, in a 
real sample not all the ferromagnetic slabs will have the 
same width determined by Ij48|l . Instead, we can con- 
sider, as in j23|, a Gaussian distribution for the widths, 
in which the mean width is given by ijisll . Figure El shows 
the behaviour of {qy) as a function of temperature for a 
Gaussian distribution whose mean standard variation is 
5% of the mean width. Note that the local minima re- 
ferring to the configurations whose number of domains 
along the x axis is large are almost suppressed, while the 
minima correspondent to a small number of domains are 
more pronounced. 

This prediction for the qualitative behaviour of {qy) is 
a consequence only of the competition between the inter- 
actions involved (which generates the semi-ellipsis at the 
phase space) and of the geometry considered (since the 
local minima appear when the modulation length tt/qo 
is "commensurate" to the slab's width d, favouring the 
modulation to lie only along the width's direction). It 
does not depend on any particular aspect of the model, 
and is expected to hold even in a vectorial model. Unfor- 
tunately, until this moment, there are no experimental 
data available to verify such prediction. 

In figure El we compare the energy difference of the 
bubble phases to the energy of the striped ones, but with- 
out the local energy minima, to make the plot easier to 
read. It is clear that the former is always greater than the 
latter, what makes one expect to not find bubbles inside 
the ferromagnetic terraces. In fact, the Magnetic Force 
Microscopy (MFM) images of MnAsiGaAs films do not 
show configurations like bubbles, but rather structures 
similar to the stripes predicted by our model and pre- 
sented in figure Eb- Moreover, as discussed in US [13, 
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Figure 10: Free energy difference between ordered and disor- 
dered states associated to the striped (full line) and bubble 
(dashed line) phases for MnAs:GaAs films as a function of 
temperature (in Celsius degrees). The lines correspondent to 
the configurations for which the y direction is non-modulated 
were removed for the sake of clearance. 



the MFM images also suggest that when the magnetiza- 
tion hes along the z direction, the ferromagnetic terraces 
are divided in 2 or 3 domains, and not in the wider range 
of 1 to 10 domains predicted by our model. Of course, 
there are several ingredients lacking in our model to make 
it more realistic, like the vector nature of the magnetiza- 
tion and the presence of topological defects that may take 
into account these features (indeed, Garel and Doniach 
showed that the occurrence of dislocations in the case of 
the infinite slab can melt the ordered phase). Nonethe- 
less, our simple model, taking into account only the na- 
ture of the interactions and the geometry involved, is 
able to describe some general qualitative features of the 
system. 



V. CONCLUSIONS 

In this work, motivated by the recent experiments that 
show the morphology of the MnAsiGaAs ferromagnetic 
stripes, we developed a general theory that describes the 
magnetic phase transitions at a dipolar Ising ferromagnet 
slab with finite thickness and finite width. We showed 
that although the modulated phase occupies a smaller 
volume in the (qx, Qy) momentum space than the one oc- 
cupied by the same phase in the case of an infinite slab, 
the transition between the ordered and the disordered 
configurations is still first-order and induced by fiuctu- 
ations (a Brazovskii type transition). This happens be- 
cause all the momentum space "shrinks" in the case of 
the finite slab, due to the discretization of the momen- 
tum component in the limited direction {x), what is a 
consequence of the boundary conditions. And what is 
fundamental to the Brazovskii transition occur is not the 
absolute volume of the modulated phase, but how large 
it is when compared to the rest of the momentum space. 
So, in the case we studied there is still enough phase 
space to the fluctuations of the order parameter induce 
a transition. 



We also showed that, for the finite slab, there is the 
possibility of first-order transitions between the modu- 
lated and the disordered phases driven not by temper- 
ature, but by variation of the slab's width. Studying 
two different modulated configurations, namely, striped 
(figure [St) and bubble (figure Cb) ones, we showed that 
the first has lower energy and higher spinodal then the 
latter, what means that bubble phases are not expected 
to be observed in such materials. This is in qualitative 
agreement with MFM images realized on MnAs:GaAs 
thin films, where the domain structures inside the fer- 
romagnetic terraces are similar to the stripes of figure 

Another effect of the finiteness of the slab was that 
modulated phases characterized by different number of 
"spread" domains (referring to the semi-ellipsis minimum 
energy projection along the qx axis ) have different ener- 
gies. Thus, fluctuations of the order parameter, together 
with flnite size effects, break the high degree of degen- 
eracy of the ground state. In the case of the infinite 
slab, there was already a break of degeneracy (referring 
to the differences in energy between striped and bubble 
phases, for instance), but in the finite case, we noticed 
that it was deeper (as some degeneracies between differ- 
ent stripes are also broken). It is clear that this effect, 
in the momentum space, is a consequence of the break 
of translational invariance in the real space. However, 
for the case of the MnAsiGaAs system studied so far, we 
saw that this break is practicable undetectable, due to 
the large film thickness (hundreds of lattice parameters) 
and temperature (that is far below the spinodal limit). 

In addition, we noted that for slabs whose widths d are 
"commensurate" with the modulation wavelength n/qo, 
the most stable configuration is the one for which there 
is no modulation along the unlimited (y) direction. This 
is refiected by the occurrence of steps in the energy pro- 
files of the system. In what concerns MnAsiGaAs films, 
we noticed that such steps generate local minima in the 
free energy profile that are responsible for changes of the 
number of "spread" domains in the x direction inside the 
ferromagnetic terraces. Although MFM images reveal 
that there can be phases, inside the terraces, with dif- 
ferent number of domains (2 and 3), they show that this 
number is not so large as the predicted by our model 
(below 10). Moreover, it was not reported yet any con- 
figuration without modulation along the y direction. In 
fact, we do not expect that this actually happens, not 
only because of the precision required (the slab's width 
must be commensurate to the modulation wavelength), 
but also for the fact that the y direction is not unlimited. 
Rather, we expect that the mean modulation along the y 
direction increases with temperature and oscillates near 
the region where the "commensurate widths" occur, as 
showed in figure El Unfortunately, as of yet there is no 
experiments that can verify this behaviour. 

Hence, this model is only the first step towards a more 
complete understanding of the complex features of the 
MnAsiGaAs phase diagram and could be applied to other 



12 



similar systems where there is a competition between or- 
ganizing and frustrating interaction plus a finite slab ge- 
ometry. An improved theory should surely comprehend 
not an Ising model, but a vectorial one, even because 
it is observed, in the MFM images, configurations in 
which the magnetization lies along the x axis. Due to 
the crystalline field of the material, no components of 
the magnetization along the y direction is expected. An- 
other important thing to take into account is the topo- 
logical defects, which can play fundamental role in two- 
dimensional systems. 



ated can be described by the magnetic scalar potential 
(j){r) that satisfies the Poisson equation [2^ : 

V^^ = -47rp, 

where p, the effective magnetic poles density, is given by: 
p = -V • M . 

Hence, the magnetostatic energy is written, in the 
Fourier space, as: 
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Appendix A: Dipolar energy 



p{k) 



fc2 



-<rk. 



(49) 



Let us now apply this formalism to our specific case, 
namely, an arbitrary magnetization of a slab with thick- 
ness D and length d. Making use of step functions, it can 
be written, in the whole space, as: 



In this appendix, we explicitly calculate the magneto- 
static energy of an arbitrary configuration of the slab. 
First, let us derive a general expression to compute mag- 
netostatic energies: given a certain magnetization M(r), 
in the absence of free currents, the magnetic field gener- 



M{x) = M{x, y)e{x)d{d - x)9{D/2 - z)e{D/2 + z)z , 

(50) 

where M{x,y) is given by (EJ, as explained before. 
Hence, a straightforward calculation yields, for the ef- 
fective magnetic poles density in the Fourier space: 



[e''=-'*(-l)" - 1] 



Substituting this expression in ijiflji . we obtain the dipolar energy: 



3 mn^ky-mn' .-ky 2 



{kx + yi-K / d) [kx — WK / d) 
1 + (-!)""'+! cos(fc:rrf) 



where the summation is to be understood as involving n and n' with same parity. Therefore, there are three integrals 
to be evaluated; the one referring to ky can be rewritten as a summation and the one referring to k^ can be calculated 
analytically: 



dk. 



sm 



{kl + kl + k^) 2 
yielding the following expression for the density of dipolar energy: 



h'^ -X Zf2 



fdip — 



E 



LyDd Dd^ V a 



\ n3 ) 2^ 



ky .n^n' 



') l + (-l)""'+icos(fc^d) 



L2 t jL2 (1.2 _ W^TT^ 



which, by a simple change of coordinates, becomes expressi on 
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